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Abstract 

Current flowing is studied in magnetic junctions consisting of a ferro- 
magnetic metal (FM) , antiferromagnetic conductor (AFM) and a nonmag- 
netic metal closing the electric circuit. The FM layer with high anisotropy 
and pinned spins of the magnetic atoms in the lattice acts as a spin injec- 
tor relative to the AFM layer. To obtain resulting magnetization in the 
AFM layer, magnetic field is applied, which may be varied to control the 
magnetization. The spin-polarized current from the FM layer creates a 
torque and makes it possible to switch the magnetization. A possibility is 
shown to lower the threshold current density by the orders of magnitude 
by means of the magnetic field. 



1 Introduction 

Antiferromagnetic layers are used conventionally in spintronics to intro- 
duce additional magnetic anisotropy and pin the magnetization vector 
direction in a magnetic junction. Such an effect is based on the unidirec- 
tional anisotropy phenomenon [I]. Recently, interest has been revived to 
studying FM-AFM junctions with special attention to the current effect 
on the interface processes J2j [3j S] - Besides, experimental works have been 
revealed [3 |Hj in which the magnetoresistive effect has been observed in 
FM-AFM structures with a point contact between the layers under high 
current densities (up to 10 9 -10 10 A/cm 2 ); the effect is due, apparently, to 
the sd exchange interaction. The similar effect is observed in FM junctions 
with broad applications [7][S]. Suppositions have been voiced that such an 
effect in FM-AFM or AFM- AFM structures will have interesting features 
because of absence of demagnetization, in particular. Finally, opinions 
have been stated that the sd exchange may lead to the spin transfer from 
the conduction electrons to the lattice, as in FM junctions, and cause 
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instability with the magnetization switching. It has been supposed that 
nonlinear effects and the current-driven electron spin polarization may 
promote such an effect P [HJ1 ITT1 IT2] . 

Thus, the possibility of the current-driven exchange switching FM- 
AFM structures remains, apparently, the most incomprehensible in theory 
and experiment. The idea of the electron spin transfer to the lattice by 
the spin-polarized current has been put forward in famous works |13l 114] , 
The idea turned out to be very fruitful, so that it is interesting to know 
whether it is applicable to the FM-AFM structures. An attempt to clear 
the matter up is the goal of the present paper. 

When we say about the spin transfer "to the lattice", it is necessary 
to explain which of two (at least!) AFM sublattices spins transfer to. The 
basic works |13l I14| , naturally, do not answer the question. A physical 
picture of the spin transfer, according to |14] , consists in precession of 
the current polarization vector p around the lattice magnetization vector 
M and gradual decrease in the precession angle with motion along the 
current because of the statistical spreading of the electron velocities. Such 
a picture looks general enough and must remain for any model of the 
magnetic sublattices, including AFM. If such a picture is applied, then 
a torque T similar to that in [T3] is to be added to the right-hand side 
of the Landau-Lifshitz (LL) equation for the AFM magnetization. Note 
that the torque T reveals in LL equation for the total magnetization M, 
rather than the each of the sublattices Mi, M2. 

A formal derivation of the expression for the T moment was proposed 
by us in Refs. [151 1161 I17| . It should have in mind, however, that the 
torque originates in a very thin layer Af ~ 1-2 nm thick at the FM-AFM 
interface; such a length coincides with the precession damping lengths of 
the M and p vectors [131 114] . So T moment exists only in that thin layer, 
not the whole range of the current interaction with the magnetic lattice. 
In accordance with Refs. [1511161 [17] . we take the T effect into account by 
means of the boundary conditions. 



2 The structure and equations of motion 
2.1 Conduction electrons in the AFM layer 

The structure in study is shown schematically in Fig. \T\ The sizes along Oy 
and Oz axes are large, so that demagnetizaion is absent in that dimensions. 
For simplicity, we assume the FM layer to be pinned. This layer injects 
spin polarized along p = Mfm (here and below the hat over the vector 
indicates a unit vector). The main processes occur in the AFM layer. Let 
us consider the equations of motion for that layer. 

The electron magnetization m(a;, t) obeys the continuity equation 

drrii , . , „ rrn — rhi , . 

-Qjr + V fc J ik + ga sd [m x M]j — — - = 0, (1) 

where g is the gyromagnetic ratio, a s d is the sd exchange constant, 

J lk = Mi^-Ut - 3k) = Mi^f {Pju - enD^L} (2) 
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Figure 1: The structure in study. The FM layer is a ferromagnetic metal with 
Mfm magnetization pinned along z axis; the AFM layer with L thickness is 
an easy-plane-type antiferromagnetic conductor with sublattice magnetization 
vectors canted by 9 angle in an applied field H parallel to the anisotropy axis 
n||z; Mi | = IM2I = Mq. The electric current is closed by a nonmagnetic 
conductor NM, so that electron flux j/e flows through all the layers. 



is the electron spin current, are the partial electric current densities in 
the spin subbands, ju — j£ + j~ is the total current density, D is the spin 
diffusion constant (it is assumed to be the same in both subbands), to, = 
mMi, m = hbuP, n is the electron density, P is the spin polarization, 
fhi = fiBtiPMi, 

3/2 1 s3/2 

p _ tp - (eg - tsd) , . 

^ _ 3/2 7 } 73/2 w 

is the equilibrium spin polarization, e s d is the spin subband energy shift, 
£f is the Fermi energy relative to the bottom of the lower subband. 

The ways of the Eq. |T} transformation and simplification were dis- 
cussed in detail in Refs. |15H16llT7] , so that we present only some necessary 
results here. 

We consider only slow enough processes due to the vector M precession 
with u) <C 1/r frequencies. A typical spin relaxation time in AFM is 
r ~ 10 -12 s. Therefore, the electrons have time to follow the precession, 
so that the time derivative may be omitted. In derivation of Eq. ^ , the 
current density may be assumed to be low, namely, j/jn 1, where 
ju = enl/r ~ 2 x 10 10 A/cm 2 . This condition is assumed to be fulfilled. 
As a result, Eq. Q is reduced (see [T51 [TBI \T7\ ) to 

d 2 m m — fa ... 

a^ + ^ = a (4) 

At the layer boundaries x — and x — L the solution must obey the 
continuity conditions for the spin currents and the differences of the Fermi 
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quasilevels (ej — e F ). Then the solution of Eq. Q takes the form (see [171 ) 

j V cos V 

Am = m — m = ubuQ^u-. , — , (5) 

JZ5 I/* + cos x 

where Qfm = O'fm — Jfm)/J ^ s the current polarization in the FM layer, 
cosx = Mfm • M(0), v = Zfm/Zafm, v* = Zfm/Znm + AZfm/Zafm, 

A = L/l <gil, Z q = — are the spin resistances of the layers, L, p q , Q q 
l-Qq 

being the spin diffusion length, the resistivity and the current polarization, 
respectively, in the q = FM, AFM, NM layers. Equation ([5]) describes 
distribution of the electron magnetization across the AFM layer. 

2.2 Magnetic sublattices in the AFM layer 

We derive the equations of motion for the Mi and M2 sublattice magne- 
tizations by means of the procedure described in Ref. [18]. The desired 
equations take the following general form: 

^= S [M,xH,] + R, (i = l, 2), (6) 

where the effective fields are determined with variational derivatives of 
the AFM energy W: 

dW 

the energy is W = J w(Mi, M2) d 3 i, the energy density is w — w ex + 

W s d + W a + WH, 



fdMy 1. , S (9L 
4- ' ~'b*J + 4 (Q - a H^ 



w ex = W-L 2 ) + I(a + a')(^) +~A* 



is the sublattice exchange interaction energy density, w„d — —a s dm ■ M 
is the sd exchange interaction energy density, 

w a = -I ((n • M) 2 + (n ■ L) 2 ) - t ( (n . M) 2 - (n • L) 2 ) 

is the anisotropy energy density (the anisotropy axis n is shown in Fig.[l|, 
wh = — H ■ M is the Zeeman energy density in the external field. 
The relaxation terms may be taken in the Gilbert form [11] 



T? - a ° 
IV — , . 

Mo 



Mi x 



dt 



(* = 1, 2), 



|Mi| = [M a | = Af . 

It is convenient to introduce total magnetization M = Mi + M2 and 
antiferromagnetism vector L = Mi — M2 instead of sublattice vectors 
Mi, M2. We assume sublattices to be equivalent. Therefore, a nonzero 
magnetization M 7^ appears only under applied field H parallel to n 
axis, as shown in Fig. [TJ The 9 angle between the sublattices and n axis 
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is determined by H field with cos6> = H/H ex , H < H ex = 28 M . The 9 
angle is related with the squares of the vector lengths: 

M 2 — L 2 = 4Mo cos 9, M 2 + L 2 = 4M 2 . (8) 

It is seen that the angle is conserved with motion if the vector lengths con- 
serve. The magnetization is assumed to be low enough, M = |M| <C Mo. 
For simplicity, we neglect demagnetization and weak relativistic non- 
collinear Dzyaloshinskii interaction [18] . 

The equations for M and L vectors take the form 



dt 

9L 



dM - g [M x H^ ft M ] + g [L x H c A ff PM ] + Rl + R 2 , (9) 



^ = ji[MxL] + j[Mx Hf ff FM ] + 9 [L x YlT] + Ri R 2 . (10) 
New effective fields appear here: 

tjFM 1/ , /s9 2 M ttFM . tt 

Hefi = -(a + a )-q-5- + H + H a + H sd , (11) 



2 V '9a; 



HAFM 1/ /n 9 L „APM /,r>\ 

eff =_ 2^ a_a ^ar T ' 



with anisotropy fields 

H* M = + /3')(n ■ M)n, H f M = ±0 - /?')(n ■ L)n. 

The relaxation terms contain vector products [M x 9M/9f], [L x dL/dt], 
and [M x <9L/<9t], [L x 9M/<9t]. By transforming these products by means 
of Eqs. @ and |TUJ| with Eqs. (1X1]) and (fT2|) taking into account, we find 
that the products are of the same order of magnitude, so that we may 

ol g ~ olg- The sd exchange effective field is 



put Ri ± R 2 = R = — 9- 

Mo 



and use the Gilbert formula with 



H sd (x, t) = - g^fj, = H eq (M)M + AH sd . (13) 

The first term in the right-hand side of Eq. (|13|l is the equilibrium con- 
tribution that directed along M vector and vanishes from Eq. @. The 
second term is due to the current. It determines with the nonequilibrium 
part Am of the electron magnetization. This term coincides with the one 
calculated in Refs. |16U17| . So we may use those results here. We have 



AH sd = Mjmh. d lS(x - 0), (14) 

(15) 



7 I l\ ^ A J V i V ~ C0S X) 

h sd = (a + a )fJ,BnQ F MA — -r- 5 — 

jd \v* + cos- 3 x) 



According to Eqs. (|ll|l and (|12|) . the processes in the AFM layer differ 
substantially from those in the FM layer. In particular, even the lengths 
of the M and L vectors do not conserve in motion. The possibility of the 
electron spin transfer to the lattice similar to that in Refs. [131 114] is to 
be discussed anew. At the same time, the electron interaction with the 
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lattice is described with H 3C j field and Eqs. (|14l) (|12[) . as in the FM layer. 
Therefore, it is necessary only, that two magnetic sublattices precess in 
synchronism as in ferromagnets. 

Such a synchronism is attained easier if the sublattices are the same. 
Let the following conditions are fulfilled: a + a' 3> \a — a'\ — > and 
f3 + p' 2> \fi — P'\ —5-0. The relaxation constant is put to be small, too 
(aa — > 0). The uniform exchange dominates in the equations because 
of the following estimates [IS]: S ~ a/a 2 ~ 10 4 » f, a ~ f0" 12 cm 2 , 
a ~ 10~ 8 cm. Besides, we have typical estimates a s d ~ 10 4 , M ~ f0 2 G. 
As a result, Eqs. ([9| and (|f 0|l are reduced to 

^= 9 [MxHr fi M ]+R, f=^xL], (16) 

These equations describe precession of the magnetization vector and 
the precession-driven forced oscillation of the antiferromagnetism vector. 
Note that the vector squares M 2 and L 2 conserve in accordance with (|16[) . 
So it follows from © that the sublattice canting angle 9 conserves also in 
motion that is substantial in our model. 



3 Spin currents and boundary conditions 



Let us derive expressions for the spin currents in the lattice. According to 
Refs. [151 1161 ITT] . these expressions follow from the equations of motion. 
However, the second of Eqs. (|16[) does not contain nonlocal or singular 
terms. Therefore, the currents may be derived from the first equation 
only. We obtain from (|16p 



go. 



M x 



d 2 M 



dx 2 



d_ 

dx 



M x 



dM 

dx 



dJi 



dx 



with a magnetization flux 



M x 



dM 

dx 



(17) 



(18) 



where a — gaM has the meaning of a magnetization diffusion constant. 

Let us consider the singular term due to sd exchange in the first of 
Eqs. (fT6]l . namely, g [M x H sd ]- By substituting (fT4]l and {15]) to it, we 
have 

d~1 

fl [MxH 3d ] = ^, (19) 



dx 



where a sd exchange flux appears 



Jsd = gh sd l [m(0) x M FM ] 6(x - 0), 



(20) 



0(x) being the Heaviside step function. 

Now let us return to the electron spin current (see Eq. (HI). Here 
k = x, and the current exists in two magnetic layers, FM and AFM. 
However, the spatial gradient of the spin polarization P is absent in the 
FM layer, so that only the first term proportional to the current density j 
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remains in Eq. (2). The spin current itself contains spins collinear to the 
FM layer magnetization, i.e., 

J FM = J(-0) = M™^P FM j. (21) 

The spin current continuity conditions at x = and x = L interfaces 
give the boundary conditions necessary to solve Eqs. (|16|) . As the lattice 
is pinned in the FM layer, we have J™ = Jm(— 0) = and J™ = 
J s d(— 0) = 0. Then we obtain for x = 

J(+0) - J(-0) + Jm(+0) + J sd (+0) = 0. (22) 

Let us project Eq. (|22p to the M direction and to plane perpendicular to 
that direction. Two equations are obtained: 

m(m-j(-o)) = J(+0), 

[m x [j(-0) X M]] = J A /(+0) + J sd (+0), (23) 

which are the desired boundary conditions. At x — L boundary, the 
continuity condition gives instead of (1221) 

J(L + 0)-J(L-0)-Jm(£-0) = 0. (24) 

After projecting to M, two conditions are obtained: 

J(L + 0) = J(L-O), 3 M (L -0) = 0. (25) 



To express (|23[) and (|25|l conditions explicitly as requirements imposed 
on the desired magnetization, we multiply the conditions by M vectorially. 
Then we obtain finally 



dM(x) 



dx 



x=0 

dM(x) 



k I M FM x 



d: 



•x 



M(0)] , (26) 
= (27) 



with a parameter 



MsQfm j 
aM ev*+ cos 2 \ 
which characterizes the torque value 



k=t ^MJ_ _ (2g) 



4 The macrospin approximation 

The problem is simplified considerably if the AFM layer is thin enough, 
so that A -C 1, L <C \J ct//3 conditions are fulfilled and the magnetization 
is almost constant in that layer. Then the following expansion is valid: 

M(x) = M(0) + M'(0)a; + ^M"(0)x 2 + ..., (29) 

where the primes mean derivatives with respect to x coordinate. Thus, 
several functions of a single variable (time t) are introduced instead of a 
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single M(a;, t) function of two variables. It means physically, that a single 
large domain ( "macrospin" ) is placed in the AFM layer [13]. Such an ap- 
proximation often corresponds to experiments and simplifies calculations. 

By differentiating (|29l) with respect to x, substituting x = L and using 
boundary conditions Q26p and 1)271) . we obtain 

M"(0) = -I^M'(O) = ~ [M FM x m] . (30) 

The only term in Eq. (|16[1 that contains the second derivative may be 
transformed with (|30|) taking into account: 

o[M(Q) x M"(0)] = fe[M(0) x [M(0) x M FM ]]. (31) 

By substituting (|3 1 p into Eq. (|16[1 . we have 

+ 5 [M(0) x H'] + ^[M(0) x [M(0) x M FM ]] - 

(it -L/ 

dM(0) 



M(0) x 



dt 



= 0, (32) 



where H' = H + H FM . 

At first sight, Eq. Q32p coincides practically with that in [T3] that 
was reproduced in many publications. There is an important difference, 
however. In Eq. @ sd exchange interaction is absent because of fulfilling 
condition (|14p . so that the equation is valid, in principle, at any distance 
from x = plane where torque T originates. At a; > Af, vectors Mfm and 
M become collinear, so the torque vanishes. However, the spin current 
created in the lattice by the torque remains. 



5 Exchange instability conditions 



Let us discuss solution of the linearized dynamical equation (|32|) for small 
harmonic fluctuations AM X , AM y ~ exp(— icot) near the stationary state 
M = z. The dispersion equation takes the form 



uj~ + 2ivu> — w = 0, 



(33) 



where the following notations are used: 



Q x Qy -(- 



ak 



(34) 



-(Q X +Qy) 



ak 
La' r , 



M ■ n 



(35) 



Q x = g((H ■ M) + H™ + 4nM), Q y = g((H ■ M) + i/ FM ). 

The general instability condition Imcj > requires fulfilling inequality 
v < 0. Put, for example, the following typical parameter values: Q x ~ 
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47vgM ~ 10 11 c -1 , n y = g(H + H^ FM ) ~ 10 9 c _1 , a' G « KT 1 < 1. Then 
^ < condition is approximately reduces to the following inequality: 



1 + ^k(^- n )<°- < M > 

By solving the inequality (|36p with respect to the current density and 
taking the definition (|28p into account, we obtain the following instability 
condition: 

-(M.n)l>X a y^(l + ±). (37) 
V / e ^bQfm V " / 

Note that we discuss only one current direction when the electrons flow 
from the FM layer towards to the AFM layer. It is seen from (|37[) that 

the instability occurs only at • < under such a current direction. 

As n vector is co-directed with Mfm, it is clear that two magnetizations 
Mfm and M should be antiparallel to create instability. 

It follows from (|37l) that the threshold current density is proportional 
to the dissipation constant and squared magnetization. The latter is de- 
termined with the field and can be decreased. For example, if M ~ 10 
G is taken, then the threshold current density decreases by two orders of 
magnitude in comparison with the typical value for the FM junctions. 

6 Conclusions 

The present work has preliminary character, however, it allows to make 
some conclusions: 

• FM-AFM structure may be canted in a magnetic field, i.e., a result- 
ing magnetization can be induced. It allows to consider an exchange 
switching by means of a current-driven torque similar to that in a 
FM junction. 

• To decrease the threshold current, the AFM layer must have minimal 
dissipation and minimal value of the field-induced magnetization. 

• The instability is possible under antiparallel relative orientation of 
the FM and AFM layers. 
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